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Abstract 

We review the analytical methods of solving the stochastic equa- 
tions for barrier-type dynamical behavior in plasma systems. The 
path-integral approach is examined as a particularly efficient method 
of determination of the statistical properties. 
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1 Introduction 

A large number of physical systems, in particular in plasma physics, are 
described in terms of a variable evolving in a deterministic velocity field 
under the effect of a random perturbation. This is described by a stochastic 
differential equation of the form 

i = V(x)+V2D£(t) (1) 
where V (x) is the velocity field and the perturbation is the white noise 

(t(t) £(0> = Ht - (2) 

This general form has been invoked in several applications. More recently, 
in a series of works devoted to the explanation of the intermittent behavior 
of the statistical characteristics of the turbulence in magnetically confined 
plasma it has been developed a formalism based on barrier crossing. Previous 
works that have discussed subcritical excitation of plasma instabilities are 
Refs. P, 0, 0, 0, 0, @, 0, 0, S- 

We offer a comparative presentation of two functional integral approaches 
to the determination of the statistical properties of the system's variable x (t) 
for the case where the space dependence of V (x) is characterized by the 
presence of three equilibrium points, V (x) = 0. We will take 

V (x) = ax - bx 3 , with a, b > (3) 
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2 The functional approach 



2.1 Overview of the functional methods 

In studying the stochastic processes the functional methods can be very useful 
and obtain systematic results otherwise less accessible to alternative methods. 
The method has been developed initially in quantum theory and it is now 
a basic instrument in condensed matter, field theory, statistical physics, etc. 
In general it is based on the formulation of the problem in terms of an action 
functional. There are two distinct advantages from this formulation: (1) the 
system's behavior appears to be determined by all classes of trajectories that 
extremize the action and their contributions are summed after appropriate 
weights are applied; (2) the method naturally includes the contributions from 
states close to the extrema, so that fluctuations can be accounted for. 

There are technical limitations to the applicability of this method. In the 
statistical problems (including barrier-type problems) it is simpler to treat 
cases with white noise, while colored noise can be treated perturbatively. In 
the latter case, the procedure is however useful since the diagrammatic series 
can be formulated systematically. 

The colored noise can be treated by extending the space of variables: the 
stochastic variable with finite correlation is generated by integration of a 
new, white noise variable. 

The one dimensional version can be developed up to final explicit re- 
sult. Since however the barrier type problem is frequently formulated in 
two-dimensions, one has to look for extrema of the action and ennumerate 
all possible trajectories. It is however known that, in these cases, the behav- 
ior of the system is dominated by a particular path, "the optimum escape 
path" , and a reasonnable approximation is to reduce the problem to a one- 
dimensional one along this system's trajectory. 

2.2 The path integral with the MSR action 

We will briefly mention the steps of constructiong the MSR action functional, 
in the Jensen path integral reformulation. We begin by choosing a particular 
realisation of the noise £ (t) . All the functions and derivatives can be discre- 
tised on a lattice of points in the time interval [— T, T] (actually one can take 
the limits to be ±oo). The solution of the equation (JTJ) is a "configuration" 
of the field x (t) which can be seen as a point in a space of functions. We 
extend the space of configurations x (t) to this space of functions, including 
all possible forms of x (t) , not necessarly solutions. In this space the solution 
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itself will be individualised by a functional Dirac 5 function. 



6(x-V(x)-V2D£(tj) 



Any functional of the system's real configuration (i.e. solutions of the equa- 
tions) can be formally expressed by taking as argument an arbitrary func- 
tional variable, multiplying by this 5 functional and integrating over the space 
of all functions. 

We will skip the discretization and the Fourier representation of the S 
functions, followed by reverting to the continuous functions. The result is 
the following functional 



D[x(t)}D[k (t)]exp 



— T 



dt ( -kx + kV (x) + V2Dk£ 



The label £ means that the functional is still defined by a choice of a particular 
realization of the noise. The generating functional is obtained by averaging 
over £. 



J dt (-kx + kV (x) + iDk 2 ^j 



Z= (Z s ) = J D [x (t)} D [k (t)} exp 
We add a formal interaction with two currents 

dt ( -kx + kV (x) + iDk 2 + J x x + J 2 k 



(4) 



Zj = J D[x (t)]D[k (t)]exp 



(5) 

in view of future use to the determination of correlations. This functional 
integral must be determined explicitely. The standard way to proceed to the 
calculation of Zj is to find the saddle point in the function space and then 
expand the action around this point to include the fluctuating trajectories. 
This requires first to solve the Euler-Lagrange equations 



k = —ak + 3bkx 2 — J\ 
x = ax — bx 3 + 2iDk + Ji 



(6) 



The simplest case should be examined first. We assume there is no determin- 
istic velocity [a = b = 0) in order to see how the purely diffusive behavior is 
obtained in this framework 

x — 2iDk = J 2 
k = -Ji 
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The equations can be trivially integrated 

x (0) (t) = I dtA n {t,t')J 1 {t')+ [ dtA 12 {t,t')J 2 {t') 

J-T J-T 

k^(t) = [ dtA 21 {t,t')J 1 {t') + [ dtA 22 (t,t')J 2 (t') 

J-T J-T 



where 



A n (t,t') 

A12 (t,f) 

A21 (t,f) 

A 22 (M') 



2iD [te (t' -t) + t'Q (t - t')} 
&(t- t') 

e (t' - 1) 




with the symmetry 



Ay (t, t') = Aji (t',t) 

The lowest approximation to the functional integral Zj is obtained form 
this saddle point solution, by calculating the action along this system's trajec- 
tory. We insert this solutions in the expression of the generating functional, 
for V (x) = 



D[x(t)]D[k (t)]exp 



dt ( — kx + iDk 2 + JiX + J 2 k 



x (0) k (0) 



exp 



\% [ dt f dt'Ji (t) Aij (t,f) Jj(t) 

^ J-T J-T 



The dispersion of the stochastic variable x (t) can be obtained by a double 
functional derivative followed by taking Jj = 0. We obtain 



(*(*)*(*')> 



5 2 



4°)5Ji(t)5Ji(t') eXP 
D min (t, t') 



dt 



dt'Ji (t) A^ (t,t') Jj (t) 



Jl,2=0 



which is the diffusion. The same mechanism will be used in the following, 
with the difference that the equations cannot be solved in explicit form due 
to the nonlinearity. 

In general the nonlinearity can be treated by perturbation expansion, if 
the amplitude can be considered small. This is an analoguous procedure as 
that used in the field theory and leads to a series of terms represented by 
Feynman diagrams. We can separate in the Lagrangian the part that can 



5 



be explicitely integrated and make a perturbative treatment for the non- 
quadratic term. This is possible when we assume a particular (polynomial) 
form of the deterministic velocity, V (x). Obviously, this term is k (t) x (t) 3 
in Eq.(jSJ. The functional integral can be written, taking account of this 
separation 



5j = exp 



T r 

/ / dt(-b) 
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5 



/n i5J 2 (t) iSJx (t) iSJi (t) iSJi (t) 

where the remaining part in the Lagrangian is quadratic 



if (7) 



if = \ D[x(t)]D[k(t)]exp 



i I dt ( — kx + akx + iDk 2 + J\x + J 2 k 



The Euler-Lagrange equations are 



k = —ak — J\ 

x = ax + 2iDk + J 2 

The solutions can be expressed as follows 



x (t) = xq exp (at) + / dt Ay (t, t) Jj {t') 

Jo 



with 

An (t,0 



A 21 (t,tf) 
A 12 [t,t') 
A 22 (t,t') 



k(t) 



iD 



dt'A 2j (t, t') Jj (t') 



exp (at) 



(9) 



(10) 



x {[exp (-2a*) - 1] 6 (f -t) + [exp (-2a*') - 1] 9 (t - t')} 
x exp (at') 

exp (—at) (*' — t) exp (at') 
exp (a*) O (t — *') exp (—at') 




The form of the generating functional derived from the quadratic part is 

if = exp ii f dtx exp (at) J x (t) + % - ! dt f dt' Ji (t) A tj (t, if) Jj (t') 
I Jo 2 J J 
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The occurence of the first term in the exponent is the price to pay for not 
making the expansion around xq. However, such an expansion would have 
produced two non-quadratic terms in the Lagrangian density: (— 3bxo) s 2 k 
and (—6) e 3 k. This would render the perturbative expansion extremly com- 
plicated since we would have to introduce two vertices : one, of order four, is 
that shown in Eq.(0) and another, of order three, related to the first of the 
nonlinearities mentioned above. 

Even in the present case, the calculation appears very tedious. We have 
to expand the vertex part of Zj as an exponential, in series of powers of 
the vertex operator. In the same time we have to expand the exponential in 
Eq. ljllj) as a formal series. Then we have to apply term by term the first series 
on the second series. The individual terms can be represented by diagrams. 
In this particular case we have a finite contribution even at the zero-loop 
order (the "tree" graph). It is however much more difficult to extract the 
statistics since we will need at least the diagrams leaving two free ends with 
currents J\. 

In the case we examine here, the perturbative treatment is not partic- 
ularly useful since the form of the potential (from which the velocity field 
is obtained) supports topologically distinct classes of saddle point solutions 
and this cannot be represented by a series expansion. 



2.3 The Onsager-Machlup functional 

To make comparison with other approaches, we take = and integrate 
over the functional variable k. 



D [x (£)] D [k (t)} exp < / dt 



Dk 2 + i -kx + kV 



D [x (t)] D [k (t)) exp<- dt 



Dk + 



ix — iV 
2VD 



t lix — iV 

dt+ 

AD 



Z=Ni I D[x(t)]exp [ dt(x-V S ^ 



AD 



In other notations 



where 



Z=N X J D[x (t)) exp 



S_ 
~D 



S 



1 /• \ 2 

dt- [x-V 



This is Eq.(25) of the reference Lehmann, Riemann and Hanggi, PRE62(2000)6282. 
In this reference it is called the Onsager-Machlup action functional and the 
analysis is based on this formula. 

However, we can go further and we will find inconsistencies. We now take 
account of the fact that the velocity is derived form a potential 



V [x (t)} 



dU [x (t)] 

dx 
-U> [x (*)] 



S = J dt~ (x + 2xU' + U' 2 ^ 



dt 



1 ( 2 t tI2 "\ ldU 

4 V + U ) + 2M 



1 



-[u(t)-u(o)\ + ^ dt 



T 



This leads to the form of the generating functional 

-K (x^ t, Xi^ tj^j 

u iui 

eX P [—2D 

with 



exp 


U(T)~\ 
2D 


exp 


£7(0)1 
2D 



K (x, t, Xj , tij 



S 



dt 



(12) 



These are almost identical to the formulas (2-5) of the reference [32J (except 
that e — > 2D). Also, it is quite close of the Eqs.(7a-7c) of the ref . pH] . 
However there is an important difference. 

There is a term missing in Eq. (|T2|) which however is present in the two 
above references. The full form of the action S, instead of Eq. (fT2*|) is 



S 



w 







dt 






jj,2 




T 





D 



U" 



This term comes from the Jacobian that is hidden in the functional 5 
integration. 
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2.4 Connection between the MSR formalism and Onsager- 
Machlup 

In our approach the most natural way of proceeding with a stochastic differen- 
tial equation is to use the MSR type reasonning in the Jensen reformulation. 
The equation is discretized in space and time and selectd with 5 functions 
in an ensemble of functions (actually in sets of arbitrary numbers at every 
point of discretization). The result is a functional integral. There is however 
a particular aspect that needs careful analysis, as mentioned in the previous 
Subsection. It is the problem of the Jacobian associated with the 5 functions. 
This problem is discussed in Ref.jH]. 

The equation they analyse is presented in most general form as 



d<t>j (t) 
dt 



-(r, 



0)jk 



5H 



+ v j [<f>(t)] + e j 



where the number of stochastic equations is N , H is functional of the fields, 
Vj is the streaming term which obeys a current-conserving type relation 



— Vj [<f>]exp{-H [<!>]} = 



The noise is 6j. 

The following generating functional can be written 



Z e = J D [fy (t)] exp J dt [ifa (t)] J] 8 + K . [0 {t)] _ 0. ) j [0] 



.i-t 



the functions lj (t) are currents, 



#. [0 ( t)] = _(r o ) 



_5H_ 

H k (t) 



+ Vj [0 (t)} 



and J [(j)} is the Jacobian associated to the Dirac 5 functions in each point of 
discretization. 

The Jacobian can be written 



J = det 



Up to a multiplicative constant 



J = exp Tr In 



d_ SK\ 5Jt- 1') 
di + j0J p {t - V) 
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or 



Since the operator (J^) 
after taking the trace 





The factor 1/2 comes from value of the G function at zero. 

In the treatment which preserves the dual function associated to <f) in 
the functional, there is a part of the action 



Then a and a of the same coupling term from <pK [(f)} close onto a 
loop. Since G-7, is retarded, all these contributions vanish except the one 
with a single propagator line. This cancels exactly, in all orders, the part 
coming from the Jacobian. 

Then it is used to ignore all such loops and together with the Jacobian. 

We can now see that in our notation this is precisely the term needed in 
the expression of the action. 



<t>K[<t>] 




(pit) -> x{t) 
K^t)] - U'[x(t)} 




and the action (fT2"j) is completed with the new term 




Now the generating functional is 



exp 



U(T) 



2D 



K {x^ t, Xi, t{) 



Z 



exp 



2D 



with 




D [x it)} exp 
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and 




(13) 



W 




(14) 



Now the two expressions are identical with those in the references cited. 
This will be the starting point of our analysis. 

In conclusion we have compared the two starting points in a functional 
approach: The one that uses dual functions x it) and k (t) , closer in spirit to 
MSR; And the approach based on Onsager-Machlup functional, traditionally 
employed for the determination of the probabilities [31] , [22] • Either we keep 
k (t) and ignore the Jacobian (the first approach) or integrate over k (t) and 
include the Jacobian. The approaches are equivalent and, as we will show 
below, lead to the same results. 

A final observation concerning the choice of one or another method: in 
the MSR method, the trajectories include the diffusion from the direct so- 
lution of the Euler-Lagrange equations. In the Onsager-Machlup method 
the paths extremizing the action are deterministic and the diffusion is intro- 
duced by integrating on a neighborhood in the space of function, around the 
deterministic motion. 

3 The transition solutions (instantons) 
3.1 Numerical trajectories 

The equations for the saddle point trajectory are in complex so we extend 
also the variable in complex space 



The obtain a system of four nonlinear ordinary differential equations which 
can be integrated numerically. 

A typical form of the solution x (t) is similar to the kink instanton (i.e. the 
tanh function). The function x (t) spends very much time in the region close 
to the equilibrium point; then it performs a fast transition to the neighbour 
equilibrium point, where it remains for the rest of the time interval. 



x 



k 
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0.1 




-0.9 - 



2 4 6 8 10 12 14 16 

time 

Figure 1: The kink-like instanton connecting a position close to the left 
(stable) equilibrium point to the middle x = (unstable) equilibrium point 
in the effective potential — W(x) 



3.2 Elliptic functions instantons (from Onsager-Machlup 
action) 

The action functional Eqs. (jl3|) and ()14j) leads to the following differential 
equation (which replaces Eqs.(jBJ)) 

l - x = (^pj x 5 + (-2ab) x 3 + - 3Db^} x 

Multiplying by x and integrating we have 

x = ± [b 2 x 6 + {-2ab) x 4 + (a 2 - 6Db) x 2 + cj 1/2 (15) 

We are interested in the functions x (t) that has the following physical prop- 
erty: they stay for very long time stuck to the equilibrium points and perform 
a fast jump between them at a certain moment of time. Then we can take 
Ci = 0. The solution can be obtained form the integration 

= = ±bt + c 2 



^ e + (-2a/b)e + a 2 /b^-6D/b 
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-0.1 1 ' ' ' ' ' ' ' ' — 1 

2 4 6 8 10 12 14 16 

time 

Figure 2: The kink-like instanton connecting the middle x = (unstable) 
equilibrium point with a position close to the right (stable) equilibrium point 
to in the effective potential —W(x) 

The upper limit x q will be specified later. For the next calculation it will be 
taken as the smallest of the roots of the polynomial under the square root. 

The details of the calculations in terms of elliptic functions can be found 
in Byrd and Friedman |42j . 

The roots of the forth degree polynomial will be noted 

e 4 + (-2o/6) e + a 2 /b 2 - 6D/b = - (£ 2 - & " & " 
where 

£i = \ A*i + a> 2 £ 3 = -y/ ai - 

such as to have £ x > £ 2 > £3 > £4; then we will use x q = £4. The notations 

are 

oi\ = cb/b , «2 = \/6D/b 



13 



The following substitutions are required 



sn 2 u = 



k 2 = 





(£4 - 







£4) (£3 - 









£4) 




£3) & - 

2 


£4) 




- £3) (£2 


-0 



5 : 

a new variable is introduced identifying the lower limit of the integral, x 

a 2 = f^| 1 >l 
si £3 



ip = amui = arcsm 



fo-O (^4-^) 



(£1 - £4) (£3 - a;) 



sn«i = smip 
The integral can be written 



M^i-o (e 2 -o (&-o &-o 



P 1 l-a 2 sn 2 (w) , 
f- / —du 



, -sn 2 (w) 

This integral can be expressed in terms of elliptic functions. We take p = ^| 
1 - a 2 sn 2 (u) 



-du 



o 1 — psri 2 (u) 

a 6 u + 3a A (p 2 - a 2 ) V, + 3a 2 (p 2 - a 2 ) 2 V 2 + (p 2 - a 2 ) 3 V 3 



p 

Here the notations are 



^1= A — sfVr= n (^ 2 ^) 



V 2 



du 



(1 — p 2 sn 2 {u)) 2 



[p 2 E (u) + (k 2 - p 2 ) u+ 



2 (p2 - 1) (jfe2 _ ^2) 

9, 9 ^9 a „, 9\ / 9,\ p 4 sn (w) cn (u) dn (u) 
+ 2p 2 A; 2 + 2p 2 -p A - 3k 2 ) IT (ip, p 2 ,k)- y 2 2 { , 
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and 



4(1 -p 2 )(k 2 -p 2 ) ^ Vo+ 
+2 (p 2 A; 2 + p 2 - 3A; 2 ) V 1 + 
+3 (p 4 - 2p 2 A; 2 - 2p 2 + 3A; 2 ) V 2 + 
p A sn (u) cn (u) dn (u) 



(1 — p 2 sn 2 (u)Y 



V = J du = u = F(ip,k) = J —= 



dip 



a/1 — k 2 sin 2 ip 



The symbol U(ip,p 2 ,k) = U(u,p 2 ) represents the Legendre 's incomplete 
elliptic integral of the third kind and ip = am (u) is the amplitude of u. The 
symbols sn, cn, dn represent the Jacobi elliptic functions. 



3.3 Typical instanton solutions 

There are several well known examples of instantons. They appear in physical 
systems whose lowest energy state is degenerate and the minima of the action 
functional (or the energy, for stationary solutions) are separated by energy 
barriers. Instantons connect these minima by performing transitions which 
are only possible in imaginary time (the theory is expressed in Euclidean 
space, with uniform positive metric). It is only by including these instantons 
that the action functional is correctly calculated and real physical quantities 
can be determined. 

From this calculation we can obtain the explicit trajectories that extrem- 
ize the action functional and in the same time reproduce the jump of the 
system between the two distant equilibrium positions. These trajectories 
will be necessary in the calculation of the functional integral. However, since 
we have eliminated the external currents and integrated over the dual func- 
tional variable k (t) we cannot derive the statistical properties of x (t) from 
a generating functional. 



3.4 Approximations of the instanton form 

In the approach based on the Onsager-Machlup action the instanton is not 
used in its explicit form (elliptic functions) in the calculation of the action. 
The reason is that the result can be proved to depend essentially on local 
properties of the potential V (x). This will be shown later. 

In the approach with dual functions, one can reduce the instanton to 
its simplest form, an instantaneous transition between two states, a jump 
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appearing at an arbitrary moment of time. Using this form as a first approx- 
imation we will calculate the solutions of the Euler-Lagrange equations and 
then the action. 



4 Fluctuations around the saddle point (in- 
stanton) solution 

4.1 The expansion of the action functional around the 
extrema 



Using the Onsager-Machlup action we have 

U(x)-U (x^ 



P (x,t;xi,ti) = exp 



2D 



K (x, t, Xj i) tj) 



The new function K has the expression 

K (x,t;Xi,ti) = J D[x (r)]exp (--^ jf dr 



x 



+ W(x) 



(16) 



The integrand at the exponent can be considered as the Lagrangean density 
for a particle of mass 1/2 moving in a potential given by 



potential 



W(x) 



In a semiclassical treatment (similar to the quantum problem, where h is 
the small diffusion coefficient D of the present problem), the most important 
contribution comes from the neighborhood of the classical trajectories, x c (r) 
that extremalizes the action S. 

The "classical" equation of motion is 



2 

X c (t 



1 .. 



dW (x) 



dx 



(17) 



To take into account the trajectories in a functional neighborhood around 
x c (r) we expand the action to second order introducing the new variables 



y(r) = x (r) - x c (r) 
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This gives 

K (x,t;xi,ti) = exp 



jjS c (x, t, Xj, tj) 



rv(t)=o (If 
/ Z%(r)]exp -- / 



x 



oxp : / (It 



\y +\y 2 w" (x c (r)) 



The deviation of the action from that obtained at the extremum x c (r), can 
be rewritten 



5S = \ [ dr y ( r ) 



y\j) 



The functional integration can be done since it is Gaussian and the result is 

1 



[det {-\£i + W"{x c {r)))] 1/2 



In order to calculate the determinant, one needs to solve the eigenvalue prob- 
lem for this operator 



Vn (t) = XnVn (t) 



withe the eigenf unctions verifying the conditions 

y n (U) = y n (t) = 



Jti 



dry n (r) y m (r) = S r , 



The formal result for K is (also Van Vleck) 

K (x,t;xi,ti) = N 



(EL 



1/2 



exp 



Sc 

D 



where N is a constant that will be calculated by normalizing P. Another way 
to calculate N is to fit this result to the known harmonic oscillator problem. 

It has been shown (Coleman) that the factor arising from the determinant 
can be written in the form 



N- 



1/2 
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where the function ijj is the solution of 



~2d^ 

with the boundary conditions 



+ W"(x c (t))V = 



i/>(U) = 



In the case where there are degenerate minima in W (x) the particle can 
travel from one minimum to another. These solutions are called instantons. 
Consider for example the potential with two degenerate maxima of —V (x) 
at ±a and with a minimum at x — 0. We want to calculate the probability 
P(a,t/2;-a,-t/2). 

The classical solution connecting the point —a to the point a is a kinklike 
instanton. The energy of this solution is exponentially small 

1 -2 

E = -x -W 
4 

~ 2W^'a 2 exp \-t(2W^) 1/2 

This solution spends quasi-infinite time in both harmonic regions around ±a 
where it has very small velocity; and travels very fast, in a short time At 
between these points (this is the time- width of the instanton). 

The special effect of the translational symmetry in time is seen in the 
presence of the parameter representing the center of the instanton. It can 
be any moment of time between —t/2 and t/2. This case must be treated 
separately and we note that this corresponds to the lowest eigenvalue in the 
spectrum, since the range of variation of the coefficient in the expansion of 
any solution in terms of eigenfunction is the inverse of the eigenvalue 



Sc r , 



The widest interval, for the variation of the center of the instanton, must be 
associated with the smallest eigenvalue and this and its eigenfunction must 
be known explicitely. Instead of a precise knowledge of the lowest eigenvalue 
and its corresponding eigenfunction we will use an approximation, exploiting 
the fact the function x c (t) is very close of what we need. 
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We start by noting that x c (t) is a solution of the eigenvalue problem for 
the operator of second order functional expansion around the instanton. This 
eigenfunction corresponds to the eigenvalue 

~^x c (t) + W"(x c (r))x c (t) = 
and has boundary conditions 

x c ~ exponential small 

very close to 0, which is be the exact boundary condition we require from 
the eigenf unctions of the operator. So the difference between x c (t) and the 
true eigenfunction are very small. Since x c (t) corresponds to eigenvalue we 
conclude that, by continuity, the true eigenfunction will have an eigenvalue 
Ao (t) very small, exponentially small. Then the range of important values of 
the coefficient Cq is very large and the Gaussian expansion is invalid since the 
departure of such a solution from the classical one (the instanton) cannot be 
considered small. 

The degeneracy in the moment of time where the center of the instanton is 
placed (i.e. the moment of transition) can be solved treating this parameter 
as a colective coordinate. The result is 



K (a, t/2; —a, -t/2) 



Ao (f) 



AnDip (t/2) 
x exp 



t/2 I f 



DS [ Xl (t - B)\ 

the parameter r in the expression of the instanton solution shows is the 
current time variable along the solution that is used to calculate the action. 
The integartion is performed on the intermediate transition moment 0. We 
also note that 

S [x/ (t - 9)] = hmS [a;/ (r)] = S 

since, except for the very small intervals (approx. the width of the instanton) 
at the begining and the end of the interval (— |, |), the value of the action 
S is not sensitive to the position of the transition moment. 



K (a, t/2; -a, -t/2) 



\ A (t) " 


1/2 ^ 


An Dip (t/2)_ 





.% V ' 2 , f s 

- ) t exp — — 



AuD J 1 \ D 



The instanton degeneracy introduces a linear time dependence of the prob- 
ability. 
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In general, for a function U (x) that has two minima separated by a barrier 
(a maximum) the potential W (x) calculated form the action will have three 
minima and these are not degenerate. The inverse of this potential, — W (x) , 
which is appears in the equation of motion, will have three maxima in general 
nondegenerate and the differences in the values of — W (x) at these maxima 
is connected with the presence of the term containing D. Since we assume 
that D is small, the non- degeneracy is also small. The previous discussion 
in which the notion of instanton was introduced and K was calculated, take 
into consideration the degenerate maxima and the instanton transition at 
equal initial and final W. 

Let us consider the general shape for — W (x) with three maxima, at 
x = b, and a. The heigths of these maximas are 

DU" 

— ^ , a = 6,0, a (18) 
U'l > 0, VI > 0, Ltf < 

This is because the extrema of — W (x) 

d [-W (x)] _ Q 
dx 

coincides according to the equation of motion to the pointes where 

x = , or x (t) = const 
and the constant cannot be taken other value but zero 

x (t) = const = 

Then, since we have approximately that x ~ V (x) (for small D) then we 
have that at these extrema of — W (x) we have V (x) — and only the second 
term in the expression of W (x) remains. This justifies Eq. ()18|) . 
It will also be assumed that U£ > U". 

4.2 The harmonic region (diffusion around the equi- 
librium points) 

4.2.1 Contribution from the trivial fixed point solutions 

We want to calculate, on a Kramers time scale, t k ~ exp (AU/D) the prob- 
ability 

P (b, t; b,0) = P (b, | b, -0 =K(b, f -;b, -0 
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We have to find the classical solution of the equation of motion connecting 
(b, — |) with (b, |). This is the trivial solution, particle sitting at b 

x c (t) = b 

We have to calculate explicitely the form of the propagator in this case 

1 



K (x,t;Xi,ti) = N 

We use the formulas given before 

N- 1 



(EL ^ 



,V2 



cxp 



Sc 

D 



1 



(n n A n ) 1/2 [AnD^(t)f 2 
where the function ip is the solution of 



(19) 



2 dr 2 

with the boundary conditions 



+ W"(x c (r))V = 



dip 
dr 




1 



(20) 



(21) 



t=u 



We use simply U' and U" for the respective functions calculated at the fixed 
point b. 
We have 



W" = 



dx 2 
1 



U" 



(U") + terms of order D 



The equation for the eigenvalues becomes 

^j"-(U") 2 iP = 

ip (t) = ai exp (\U"\t) + a 2 exp (- \U"\ t) 
and it results form the boundary conditions 

1 



ai 
a 2 



2 |c/^ exp[ - |c/,/| ^ 
1 



2\U 



^exp[|[/"|^ 
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Then 



2\U 



- {exp [\U"\ (t - U)] - exp [- \U"\ (t - U)}} 



NOTE. This is the sinh which is obtained in the calculation of the ground 
level splitting by quantum tunneling for a particle in two-well potential. 
Now we can calculate 



1 



1 



[2\U 



1/2 exp[-±|l/"|(*-*i)] 



[4nD4> (t)] 1/2 (4ttD) 1/2 (l-exp[-2|f/"|(t-t,)]) 1/2 
It remains to calculate the action for this trivial trajectory x c (t) = b 



- f 
-f 

Jti 

- f 



dt 



1.2 

-x c + W (x c ) 



dtW (x c ) 



dt 



\U' 2 (x c ) - ^U" (x c 



Since the position of the extremum of W is very close (to order D) of the 
position where V (x) is zero, and since V = —U', we can take with good 
approximation the first term in the integrand zero. Then 

S c = -^£dtU"(x c ) = ~U"(t-t i ) 

We have to put together the two factors of the propagator and take into 
account that at b we have U" > 



K (x^ t, Xj, tj) 



1 



Sc 

D 



exp 

[4nD4> (t)] 1/2 

m 1/2 ««pH^ exp ( V (t-t.) 



\2nD 



;i-exp [-2U"(t-U)}Y 



( u" \ 

K(x,t;xi,ti) = ( J 



1/2 



(1 - exp [-217" (t - U)}) 1/2 



The contribution to the action is 



2ttD; [l-exp(-2[/^)] 1/2 

/ fj// \ 1/2 
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It should be noticed that no other classical solution exists since there are 
no turning points permitting the solution to come back to b. 



4.2.2 Contribution from the trivial fixed point solution in the dual 
function approach 

We apply the procedure described for the purely diffusive case to the case 
V (x) 7^ and for this we need the solution of the Euler-Lagrange equations. 
An approximation is possible if the diffusion coefficient is small. In this case 
the diffusion will take place around the equilibrium positions —^ajb and 
+ y/a/b. Taking the equilibrium x (t) = xq = ±^/a/b in the equation for 
k (t) we have 

k(t) = [ dt'A 21 (t,t')J 1 (t') (22) 
Jo 

where 

A 2 i (t, t') = exp [(-a + 3bxl) t] 9 if - t) exp [— (—a + 3bx 2 Q ) t'] (23) 

The symbol G stands for the Heaviside function. In the equation for x (t) we 
expand around the equilibrium position 

x = x Q + e (t) (24) 

and solve the equation 

e= (a- 3bxl) e + 2iDk + J 2 (25) 
taking account of Eq. (|2^|l 

cT 



x{t)=x + [ dt' An (t, t') Ji (If) + [ dt'A l2 (t, t') J 2 (if) (26) 
Jo Jo 



where 



A n (t, f) = (2iD) exp [(a - 3bx 2 Q ) t] (27) 
1 (exp [2 (-a + 3bx 2 ) t] - l) 6 (t' - t) 



x 



2 (-a + 3bxl) 

(exp [2 {-a + 3bx 2 ) i>] - l) 9 (t - i>) } 



2(-a + 

x exp [(a — 3bxl) t'~\ 



and 



Aia (t,t') = exp [(a - 3bx 2 ) t]G(t- if) exp [(-a + 3bx 2 ) t'] (28) 
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We note that at the limit where no potential would be present, (—a + 36xq) 

0, the propagators A become 

A n (t, t') -> 2iD [t& (f -t) + t'Q (t - t')] 

a 12 ( t ,t')^e(t-t') 

1. e. the propagators of a purely diffusive process (see [H5j). 

Using the solutions Eqs. (|22|) and (|26|) we can calculate the action along 
this path. 



5*, 



dt 



-kx + kV (x) + iDk 2 + J x x + J 2 k 



We will insert the expansion Eq. <|24jl . perform an integration by parts over 
the first term and take into account the equations, i. e. the first line of Eq. (jEJ) 
and Eq.| 



dt { --ke + kV (x ) e + iDk 2 + J^ + J x e 



+h[-V> (x )- Jx]e+ J 2 ^| 



Sj = Xq / dt Ji + 



dt (eJi + kJ 2 ) 



Now we express the two solutions, for k (t) and e (t) in terms of the propa- 
gators A 



Sj = x [ dtJ x (t) + ~ / dt / dt'Ji (t) A - (t, J 3 - (t') 



(29) 



where summation over i, j = 1,2 is assumed, and A 2 2 (t,t') = 0. 

We now dispose of the generating functional of the system when this is 
in a region around xq = ±y/a/b, the fixed equilibrium points. To see what 
is the effect of the diffusion in this case we calculate for the variable x (t) the 
average and the dispersion. 



(x(t)) 



i8J x (t) 



Xq 



Jl,2=0 



which was to be expected. And 
(x(t)x(t')} = 



iSh (t) iSJi (?) 
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Jl,2=0 



+ 



X 



iS Ji (t) iS Ji (t) " 
—An (t,l/)exp(iSj) 



+ 1 / dt'An (t, f ) Ji (0 + I / di'Ai 2 (t, J 2 (0 + \f dt"J 2 (t") A 21 (t", t) 
^ Jo ^ Jo ^ Jo 

+ ± jf dt"A u (f , t") J, (t") + ±£ dt"A 12 (f , t") J 2 (t") + ±£ dt.J 2 (t) A 21 (t, t>) 



x 

x exp (iSj) 

This gives the result 

(x(t)x(t')) = x 2 Q + exp [(a — 3bxl) t] D (30) 

1 



x 



2 {-a + 3bx 2 ) 
1 



(exp [2 (-a + 3bx 2 ) t] - l) (t' - i) 



(exp [2 (-a + 3bxl) t'] - l) (t - t') 



2 (-a + 3bx 2 
x exp [(a - 3bxl) f] 

In the absence of the potential a — > and 5^0, this is simply 

(x (t) x (t')) = Xq + -D min (t, *') 

i.e. the diffusion around the position x . In the present case, we note that a — 
3bx 2 Q = V (x ) = -U" (x ) with U' ' = U" (x ) > 0. Fixing the parameters, 
we take t > t' and obtain 

(x (t) x (If)) =4 + D exp (-U»t) -j^ [exp (2£tf f ) - 1] exp (-0* f ) 
For t — t' (the dispersion) we obtain (with — > 1/2) 

(t?) = 4+ (J^j [1 - exp (-2Itf f)l (31) 



It is straightforward to calculate the higher order statistics for this pro- 
cess, since the functional derivatives can easily be done. We have to remem- 
ber that this derivation was based on the approximation consisting in taking 
the equilibrium position in the Euler-Lagrange equation for k (t). 
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Figure 3: The probability distribution in the harmonic region around the left 
equilibrium stable point. 




time 



Figure 4: Comparison between the analytical formula Eq. (|31J) (continuous 
line) and numerical integration of the stochastic equation. 
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4.3 Assembling the instanton-anti-instanton solutions 

4.3.1 The formulation based on Onsager-Machlup action 

However, since W (0) — W (b) and W (a) — W (b) are small quantities and 
we can suppose that the instantons, even if they are not exact solutions of 
the equations of motion, can give a contribution to the action. The instan- 
tons connects the points of maximum not of —W (x) (because they are not 
equal) but of a different potential, a corrected —W (x) to order one in D 
that has degenerate maxima. This potential will be called (x) and it 
will be considered in the calculation of the contribution of the instantons and 
antiinstantons. 

W {0) (x) = W (x) - 5W (x) 
5W (x) = , x < x m 

= W - W b , x > x m 

where x m is the point corresponding to the minimum situated between the 
two maximas. 

It is introduced the family of trajectories xi (r — to) which leaves x = b 
at time — oo and reach x = at time oo. They have all the same energy 

E = -W b 

and the classical action <S& . 

We now consider the travel from x = b to x = made by an instanton 
xj (r — to) with the center located at time t ] next the return made by an 
antiinstanton which is actually an instanton xj (ti — r) starting from x = 
and going to x = b with the center located at time to- With these two 
instantons we create a single classical solution 

, ,_fx I (r-to) r<±(t + ti) 

The contribution to the action of this assambled solution is 

SjA^toih) ~ (Wo-W b )(t 1 -to) + W b t + 2S bo 

r. 1 2 

xi (r - * )J 



-2 / dr 

J(t 1 +t )/2 2 

The contributions in this formula comes from the potential energy and the 
kinetic energy along the trajectory. We should remember that the action is 
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the integral on time of the density of Lagrangian, where there is the kinetic 
energy term and minus the potential energy. 



J dt^-x 



If the particle would have remained in x = b imobile for all time (— 1/2, t/2) 
then the contribution from the potential would have been 

W(b)t 

The instanton spends (ti — t ) time in the point before returning to b. 
Then it accumulates the action equal with the difference in potential between 
b and multiplied with this time interval. 

[(-W b ) - (-W )} (*! - t ) = (W - W b ) (h - t ) 

(since (-W b ) > (-Wo)). 

Define the kinetic energy and the energy 



then 



f* 1 .2 

K = -x (t)dt 
J t 2 

E = l -x 2 (t) + [-W (x)} 
x(t) = y/2(E + W (x)) 



K = f (E + W(x))dt= f E + W ^ dx 

J t Q J t X (t) 

= I f y/2(E + W (x))dx 

^ Jxo 



In this formula we have to replace the expression of the trajectory xia (t) 
and integrate. 

In our case the energy has the value of the initial position. Here the 
velocity is zero and the potential is — W (b) 

E\ x=b = -W(b) 

and the potential is actually the current value of the zeroth -order potential 

E + W (x) = [-W (b)} - [-W {0) (x)] 
= W {0) (x)-W(b) 
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If we want to calculate the integral of the kinetic energy along the trajec- 
tory, we have to consider separately the intervals where the kinetic energy 
is strongly determined by the velocity, i.e. the region where the instanton 
transition occurs, from the rest of the trajectory where, the velocity being 
practically zero, the potential is a better description and can be easily ap- 
proximated. The approximation will relay on the fact that the particle is 
practically imobile in b or in 0, after the transition has been made. So we 
will use both expressions for the kinetic energy 



*i+£o o 

K\ x x z° b = / 2 dt-x\t) + / dx[W^(x)-W b ] 
Jo * Jb 



1/2 



The terms must be considered twice and the interval of integration can be 
extended for the region of transition, since in any case it is very small 



o 

2 I dx [W (0) (x) - W b ] 



1/2 - 2 r dr\x (r - t ) 



The term S&q is 



S b0 = f dx[wW(x)-W b ] 
Jb 



1/2 



Using an approximation for the form of the instanton, it results 

SjA^t^h) = (W -W b )(t l -t ) + W b t + 2S b0 
*S,OT) 1/2 



where 



and 



cxp 



d 2 W 



(2W^) 1/2 (2A m0 -t 1+ t )] 



dx 2 



x=0 



*m0 



dx 



[WW (x) - W b ] 



1/2 



(x) - W b 



1/2 



The local harmonic approximation to is 



(0) 



W, 



(0) 



W b + ^fx 2 x > x 
W b + ^f(x-b) 2 x<x 
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NOTE Since the trajectory xja (r;io,*i) is n °t an e:ract solution of the 
equation of motion in the potential W, the expansion of the action S will 
not be limitted to the zeroth and the second order terms. It will also contain 
a firts order term, (5S/5x)\ x=Xc ^ . It can be shown that this contribution is 
negligible in the order O (D). 

The calculation of the contribution to the functional integral from the sec- 
ond order expansion around xja is done as usual by finding the eigenvalues of 
the determinant of the corresponding operator. As before, the product of the 
eigenvalues should not include the first eigenvalue since this is connected with 
the translational symmetry of the instanton solution. This time there will 
be two eigenvalues, one for to (the transition performed by the 6 — > instan- 
ton) and the second for t\ (the transition prtformed by the antiinstanton, or 
the transition — > b). Another way of expressing this invariance to the two 
time translations is to say that the pair of instantons has not a determined 
central moment and, in addition, there is an internal degree of freedom of 
the breathing solution, which actually is this pair instanton-antiinstanton. 

To take into account these modes, whose eigenvalues are zero, we need to 
integrate in the functional integral, over the two times, to an d t±. 

The measure of integration for the two translational symmetries to an d 
ti is 



This quantity is the Jacobian of the change of variables in the functional 
integration over the fluctuations around the instanton solution. The fluctua- 
tion that corresponds to the lowest (almost zero) eigenvalue is replaced in the 
measure of integration with the differential of the time variable representing 
the moment of transition. Then it results this Jacobian. 

It can be shown that in the approximation given by exponentially small 
terms, the contribution to the path integral of the small fluctuations around 
the classical instanton-antiinstanton xja solution is the product of the fluc- 
tuation terms around the instanton and antiinstanton separately. 



2itD 





1/2 



1/2 



2nD 



I 1/2 



Sia (t; t , ti) 



x 



cxp 



«') 



D 



30 



where {Xo) T A represent the lowest eigenvalue of the operator arising from the 
second order expansion of the action, defined on the time intervals 

x c (t) = x 7 (r - t ) , -- < r < — - — 

x c (r) = xj (ti - r) , — — < r < - 
The notation ipj and respectively ip A represents 

ipj value at ^° 1 of the eigenfunction ip that starts at t { ■ 



t 
2 



■0^4 value at - of the eigenfunction ijj that starts at ti = ^ 1 
This corresponds to the formula of Coleman which replaces the infinite prod- 

1/2 



uct of eigenvalues with 



A 



4-irDip 



where ip is calculated at the end of the 



interval of time, where ip verifies the boundary conditions (}2"T]) . 
The result is 



V2 , 



' 2' ' 



(32) 



2nD 



{2Wg) 
x exp 

t/2 



/A 1/2 



&)|a: m | WtfW, 



2tt 2 D 2 



exp 



fc 6 



^ + A 6m (2<) 1/2 + A &m (2<) 1/2 



yt/2 

x / dto / c?ti exp 

-t/2 Jt 

where 



5 



(0) Jb) 



(*! " to) + ^ exp [- (2^') 1/2 (*i - to) 



//\ 1/2 



is the lowest eigenvalue of the Schrodinger equation associated with the 
Fokker-Planck diffusion equation in the local harmonic approximation of the 
potential W in the well i. 



C — x r 



exp 



2 (2W£ 



//\l/2 



A 



m0 



(33) 



and 



(ar) - W b 



1/2 



W^ (0) (x) - W h 



1/2 
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There is a very important problem with this formula: the coefficient C 
is positive and the contribution of this part in the t 1 integration comes from 
time intervals 



f 1 - 1 < 



1 



(2Wg) 1/2 



At 



Since this is a very small time interval, it results that the contributions are 
due to states where the instanton and the antiinstanton are very close one 
of the other, which is unphysical. It will be necessary to calculate in a 
particular way this part, introducing acontour in the complex t\ — t plane, 
with an excursion on the imaginary axis. 

The calculation of the propagator (b, |; b, — |). Here (1) means 
that only one pair of instanton and anti-instanton is considered. 

We particularize the formula above using the expression for W (x) given 
in terms of U (x) . Now we have 



W 

w b 
w 



DU' b > 



DU'> 



up 



TTl/2 

u 



The expression of Sb is 

U -U b D 
S b0 = ^^ + -\n 



U' b ' (x m - b) 



Uq x m 



+ y (U^A hm + ^'A m0 ) 



u' b ' 


U" 


,->-» 2 




U" 




2 


\2ttD 



3/2 



I(t) 



x exp 



Uo-U b 
D 



+ 2A 



mO 



\TJ" 



where 



,t/2 , 

I (t) — / dto / dz exp 

J-t/2 JC 



\U»\z + ^exp(-\Ug\z) 
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4.3.2 The formulation using dual functions 



We will use the functional approach in the setting that has been developed 
by us to the calculation of the probability of transition from one minimum 
to the same minimum with an intermediate stay at the unstable maximum 
point (symmetric potential). 
The initial equation is 



We have to solve these equations, and replace the solutions x (t) and k (t) in 
the action functional. Then the functional derivatives to the external current 
J\ (t) will give us the correlations for the stochastic variable x(t). We can 
also calculate the probability that the particle, starting from one point at a 
certain time will be found at another point at other time. This will be done 
below. 

The first step is to obtain an analytical solution to the Euler-Lagrange 
equation. The method to solve these equation is essentially a successive 
approximation, as we have done above, for the diffusive motion around a 
stable position in the potential (harmonic region). We know that the classical 
trajectory must be of the type of a transition between the initial point, taken 
here as the left minimum of the potential and the final point, the unstable 
maximum of the potential. From there, the particle will return to the left 
minimum by an inverse transition. We have to calculate simultaneously x (t) 
and k (t) , but we have sufficient information to find an approximation for 
x (t) , by neglecting the effect of diffusion (the last term in the differential 
equation for x (t)). 



x(t) = v(x)+at) 



The action functional with external current added is 




dtL 



L= -k (t) x{t)+k (t) V (x) + iDk 2 (t) + J x (t) x (t) 
The equations of motion are 

x-V(x)- HDk = 




x = ax — bx 3 



This gives 



x(t) = ± 




&{l + exp[-2a(r-t )]} 



1 



1/2 
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This solution is a transition between either of the minima ±ya/6 and 0. 
It shows the same characteristics as found numerically or by integrating the 
elliptic form of the equation, in the case of Onasger-Machlup action. The 
particle spends long time in the initial and final points and makes a fast 
transition between them at an arbitrary time to- The width of transition is 
small compared to the rest of quasi-imobile stays in the two points, especially 
if a is large. Then we will make an approximation, taking the solution as 



where -ya/6 and are the initial and final positions, the indice 77 means that 
this is the calssical solution (extremum of the action) and (1) means the first 
part of the full trajectory, which will also include the inverse transition, from 
to — \/ajb. The following structure of the total trajectory is examined: 
The total time interval is between — T and T (later the parameter T will 
be identified with —t/2 for comparison with the results from the literature). 
The current time variable is r and in the present notations, t is any moment 
of time in the interval (— T, T). At time to the particle makes a jump to the 
position x = and remains there until r = t\. At t\ it performs a jump to 
the position — ^a/b, where it remains for the rest of time, until T. 




With this approximative solution of the Euler-Lagrange equations (since 
we have neglected the term with D in the equation of x (r)) we return to the 
equation for k (r). We first calculate 



We need the integration of this quantity in the inverse direction starting from 
the end of the motion toward the initial time 

[ dr (^] = [-a (T - t)} 9 (t - t ) + [-aT + 3at - 2at] (t - t) 
Jt \dxJ x w {T) 

We now introduce the second part of the motion: at time t — t\ the 
particle makes the inverse transition 





dV 



(-2a) 9 (t - r) + a0 (r - t ) 




with the similar quantities. 
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After explaining the steps of the calculation, we change to work with the 
full process, assembling the two transitions and the static parts into a single 
trajectory 

z„M = (r) + (r) 

dV_ 

dx 

and 



= (-2a) (t - r) + a© (h - r) (r - t ) + (-2a) (r - h) 

x v (t) 



Jt \dxJ x m {T) 



We have introduce the notation 

W{t) = (2a) (T — t)Q(t — ti) (34) 
+ (2aT - 3ati + at) (t - t ) (h - t) 
+ (2aT - 3a*i + 3at - 2at) (t - *) 

Then the solution of the equation for the dual variable is 

k v (r) = exp [-W (r)] \k T + j dt' [-J x (t')] exp (*')] 

According to the procedure explained before we will need to express the 
solutions as bilinear combinations of currents, so we identify 

A 2 i (t, *') = (f - t) exp [-W (t)} exp [W (t')\ (35) 
and the solution can be rewritten 

k v (r) = k T exp [-W (r)] + J dt'A 21 (t, t') J x (*') (36) 

Using these first approximations for the extremizing path (x v (r) , k v (r)) 
we return to the Euler Lagrange equations and expend the variable x as 

x (r) = x v (r) + 5x (r) (37) 

whose equation is 

5x (r) = 2iDk v (r) + (^) 5x (r) 



\dx J 



x v (t) 
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and the solution 

Sx (t) = exp 







\ / x„(t) 



B + J dt'2iDk v (t')exp 



Here B is a constant to be determined by the condition that 5x vanishes at 
the final point. We note that here all integrations are performed forward in 
time. We also notice that this expression will contain the current J\ and we 
will introduce the propagator An. Performing the detailed calculations we 
obtain 

Sx(t) = B exp[W(t)]+exp[W(t)} j dt' (2iDk T )<S) (t - if) exp[-W 

+ j dtf&u(t,if)Ji(if) 

The first two terms depend on constants of integrations and the propagator 
is 

A n (t, t') = exp [W (t)} [ dt"2iDQ (t - if') A 21 (t", if) exp [-W {t")\ (39) 

J-T 

Using the solutions (|3*7jl and (|3l)jl we can calculate the action along this 
trajectory. 



S 



dt -kx + kV + iDk 2 + J x x 



dt 



■T 

T fl ■ I 

dt —xk + —kV + J\x 



-kx + kV + 2iDk 2 + - -kx + kV + J x x 



dt 



1 

r 



Ji + k— + -kV + Jxx 



dt { -J\X + -k 



dx J 

ax 



[It becomes obvious that when the potential is linear (which means that 
the Lagrangean is quadratic the two functional variables x (t) and k (t)) the 
potential does not contribute to the action along the extremal path]. 
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Using the solutions we have 

s = 1 

dt- 



dt-Jx (t) [x v (t) + 5x (t)} 
2 



+i 



V{x v ) 



dV\ 
dx J 



X 7 



dt'A 21 (MVi(0 



It will become clear later that the parts containing constants are not sig- 
nificative for the final answer, the probability. We will focus on the terms 
containing explicitely the current J\ since the statistical properties are de- 
termined by functional derivatives to this parameter. The action is 

r-T 



S = 



dt 



+ 2^i (0 



-Ji(*K(f) 



dt' An (M'W) 



+ - 



V(x v (t))- 



dt'A 21 (*,0^i(0 



The generating functional of the correlations (at any order) is 

Z[J 1 ]=exp{S[J 1 ]} (40) 

and we can calculate any quantity by simply parforming functional deriva- 
tives and finally taking J\ = 0. 

Instead of that and in order to validate our procedure, we will calculate 
the probability for the process: a particle in the initial position x = — a/ a/b 
at time — T can be found at the final position x = x c at the time t = t c . 
(Later we will particularize to x c = —^Ja/b and t c — T). The calculation of 
this probability P can be done in the functional approach in the following 
way. 

We have defined the statistical ensemble of possible particle trajectories 
starting at x — —y/a/b at t — — T and reaching an arbitrary point x at time 
t. In the Martin-Siggia-Rose- Jensen approach it is derived the generating 
functional as a functional integration over this statistical ensemble of a weight 
measure expressed as the exponential of the classical action. Now we restrict 
the statistical ensemble by imposing that the particle is found at time t — t c 
in the point x = x c . By integration this will give us precisely the probability 
required. 

The condition can be introduced in the functional integration by a Dirac 
5 function which modifies the functional measure 



D[x (r)]-> D[x(r)]S[x (t c )-x c ] 
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and we use the Fourier transform of the 5 function 

Fa, 



P I x Cl t c , 



-T 



D [x (r)] 6 [x (t c ) - x c \ exp (S) 
— J dX exp (— i\x c ) J D[x (r)] exp [S + i\x (t c )] 



It is convenient to write 



S + iXx (t c ) = % \ drL [x (t) , k (t)] + i / drXx (r) 5 (r - t c ) 



-T 
T 



—T 



[ dr[L + Xx (r) 5 (r - t c )\ 

J-T 



We can see that the new term plays the same role as the external current 
Ji (r) and this sugests to return to the Eq. (j40|) and to perform the modifi- 
cation 

Ji(t)-> J[{t) = Ji(r) + A5 (r-t c ) 

obtaining 

p(x c ,t c ;-^,-Tj = 7^ J dXexp(-tXx c )tf x) [J x ] 



Ji=0 



Using the previous results we have 



[Ji] 



Ji=0 



exp y-Xx v (t c ) 

+i-X 2 A u (t c ,t c ] 

dV 
dx 



(41) 



V 



x 



XA21 (t c , t c 



t=t c 



Ji=0 



This expression will have to be calculated at t c = T and x c 



' 'a/b. 



V 

dV 

dx 



so the last term is zero. 



[Ji] 



Ji=0 



exp 








(42) 



Ji=0 



38 



We need the propagator A 2 i in the expression for An, 

A 21 (*', T) = 6 (T - exp [-W (If)] exp [W (T)] 
W (T) = 

Then the intration giving the propagator An is 

J dt'2iDQ (T - f) 6 (T - if) exp [-W (t')) exp [-W 
= 2iD J dt'exp[-2W (t')} 
= 2iD | exp [— 4aT + 6a£i — 6ato] J dtexp (4at) 

+ exp [— 4aT + 6ati] / dt exp (— 2ai) + 

•/to 

+ exp [— 4aT] J dt exp (4at) | 
We introduce the notation 

y(t ,ti) = 1 + 

+3 exp [-4a (T - t x )\ {exp [2a - t )] - 1 
— exp [— 8aT] exp [6a (ti — to)] 



and we obtain the expression 



iD, 



A n (T,T) = —Y(t ,t 1 ] 
2a 



Returning to 

[J x = 0] = exp 

and the probability is 



A la A 2 £> ' 

-«-4/ t- : F 



2 V 6 4a 



P (to, *i 



27 



1 

27 



dAexp 
dX exp 



zA 



2 <A V6 



[Ji = 0] 



a AD 



4a 



Y(to,h 
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We have made manifest the dependence of the probability on the arbitrary 
times of jump t ,ti. Integrating on A 



P(t ,t!) = 



ixDY 



exp 



a 2 \ 
AbDY ) 



This expression contains an arbitrary parameter, which is the duration 
of stay in the intermediate position, i.e. at the unstable extremum of the 
potantial, x — 0. Integrating over this duration, t — t\ , adimensionalized 
with the unit 1/a, for all possible values between and all the time interval, 
2T, will give the probability. In doing this time integration we will make a 
simplification to consider that T is a large quantity, such that the exponential 
terms in Y will vanish. Then Y reduces to 1 and the integration is trivial 



p2T 

P = d[a(*i-* )]P(to,*i) 
Jo 



2Ta\ i exp 

ixD 



— ) 

AbD 



This is our answer. To compare with the result of Caroli et al., we first adopt 
their notation T = t/2 and we write the potential as 



V (x) = -U' 



x) 



with 



We note that 



U(x) = --ax 2 + -bx A 
2 4 



U (x = 0) - U - 




+ 
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which allows to write 



« 3/2 

P = t , exp 



'tiD 



D 



which corresponds to the result of Caroli et al. for a single instanton-anti- 
instanton pair, where 



a 3 / 2 ~ \U, 



//| 1/2 



U' 1 



'a/b 



The only problem is the need to consider the reverse of sign due to the 
virtual reflection at the position x — 0. In the tratment of Caroli, it is 
included by an integration in complex time plane, after the second jump, t±. 
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4.4 The trajectory reflection and the integration along 
the contour in the complex time plane 



The problem arises at the examination of the Eq.()32|) where it is noticed that 
the constant C multiplying the exponential term in the exponential in the 
integrations over t and t\ is positive (see Eq. (J3H)) ). Then the dependence of 
the integrand on t\ — to makes that the most important contribution to the 
integration over the variable t\ to come from 



{2W^ 1/2 



Here the integrand is of the order 



C 

0X1 } 1 D 

This would imply, in physical terms, that the most important configurations 
contributing to the path integral would be pairs of instanton anti-instanton 
very closely separated, i.e. transitions with very short time of stay in the 
point 0. This is not correct since we expect that in large time regime the 
trajectories assembled from instantons and anti-instantons with arbitrary 
large separations (duration of residence in before returning to b) should 
contribute to K. 

This problem can be solved by extending the integration over the variable 
representing the time separation between the transitions 

z — t\ — t 

in complex. 

NOTE. The quantity 

{2W^f 2 = Lu 

is the frequency of the harmonic oscillations in the parabolic approximation 
of the potential. 



The contour C is defined as 

Im z 
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1 77" I 



t 

sC Re z ^ - - t 

The form of the trajectory xja_ connecting the point b to itself, b, assembled 
from one instanton xj and one anti-instanton xa can now be described as 
follows. 
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The trajectory starts in the point b and follows the usual form of the 
instanton solution, which means that for times 

~<r<to 

it practically conicides with the fixed position b. At to it makes the transition 

to the point 0, the duration of the transition (width of the instanton) being 

At. Then it remains almost fixed at the point 0. The problem is that there 

is no turning point at x = since there the potential — W (x) is smaller than 

at x — b. Then we do not have the classical instanton which takes an infinite 

time to reach exactly x — 0, etc. Also there is no exact instanton solution 

which starts from x = to go back to x — b. In the assambled solution 

xj + xa, there is a discontinuity at x = in the derivative, the velocity 

simply is forced (assumed) to change from almost zero positive (directed 

toward x — 0) to almost zero negative (directed toward x — b). What we 

do is to introduce a turning point, which limits effectively the transition 

solution, but for this we have to go to complex time, to compensate for the 

negative value of the potential. Making time imaginary changes the sign of 

the velocity squared (kinetic energy) and correspondingly can be seen as a 

change of sign of the potential which now becomes positive and presents a 

wall from where the solution can reflect. This is actually a local Euclidean- 

ization of the theory, around the point x — 0. 

At time (to + ^i) /2 the trajectory starts moving along the imaginary z 

axis. It reaches an imaginary turning point at t^t and returns to the point 

I o I 

Rez = xia (^), Imz = 0. 

This trajectory can be considered the limiting form of the trajectory of 
a fictitious particle that traverses a part of potential (forbidden classically) 
when this part goes to zero. 



Then the expression is 



( 






1/2 (W|) 1/2 


\2ttD) 


2tt 



U -U b 



D 

where the minus sign is due to the reflection at the turning point. 



4.5 Summation of terms from multiple transitions 

Calculation of the transitions implying the points b and a 

This is the term 




42 



and corresponds to the transitions 



Again we have to introduce a new potential replacing W. The real potential 
W (x) as it results from the Langevin equation, depende on D which is sup- 
posed to be small. Then the "expansion" in this small parameter D has as 
zerth order the new potential (x) that will be used in the computations 
below. This is 

W {0) = W-8W 



5W (x) 




where the new point x p is in the region of the right-hand minimum of — W (x). 
This point will correspond to the center of the instanton connecting to a. 
We introduce the new family of instanton solution 

xi (t - h) 

in the potential — that leave x = at — oo and reach x = a at time oo. 

With this instanton and the previously defines ones, we assamble a func- 
tion that provide approximately the transitions conncetin b to b via a. This, 
obviously, is not a classical solution of the equation of motion, although it is 
very close to a solution. 

t < \ (t + t x ) 

/_.-/ f / f ^ _ J •' ' V 'M | (*0 + h) < T < \(ti + t 2 ) 

x iiaa (tMiMM)-{ l( tl + t 2 )<r<l(t 2 + t 3 ) 

r > ^ 

These trajectories are rather artificially, being assambled form pieces of true 
solutions. They have singularities: in the middel, at (ti + t 2 ) /2 the function 
is continuous but the derivative is discontinuous; in the points where one 
instanton arrives and another must continue (i.e. (t + ti) /2 and (t 2 + t 3 ) /2, 
both the function and the derivatives are discontinuous. 

this is corrected by extending the trajectories in complex time plane. One 
defines the variables 



xi (r - 


-t ) 


xi (r ■ 


-h) 


Xi (t 2 


-r) 


Xl (t 3 


-r) 





= h 


- to 


z 2 


= t 2 


-h 


z 3 


= h 


-h 
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and the three integration contours 
C\ Im z\ 



\u"\ 



^ Rezi ^ \ -t 



C 2 lmz 2 = ^ 7T < Rez 2 < Re (f - h) 



C 3 Im z 3 



\u"\ 

o 



^ Rez 3 ^ | -t 2 



Together these contours must assure the continuity of the function xja in the 
region 0, from which the condition arises 



Xi 



Re 



t\ — to 



-Xi 



Re 



to — ti 



which gives for large separations ti — t 

x p exp (\Uq\ A 0p ) = \x m \ exp (\Uq\ A m0 ) 
The exprsssion of K at this moment is 



K, 



(2) 
ba 



t 



where 



U>> ) 4 mp \2::D 



\U, 



in \ 3 



x exp 

"t/2 



Wp -Ug-U b 

D 



+ 2|f/ 7 |(A m0 + A 0p ) 



x 



/ dt [ dz x \ dz 2 [ dz 3 exp[F( Zl ) + G(z 2 ) + F (z 3 )} 

J-t/2 JCx JC 2 JCz 



F{z) = ^exp(-K'|z) 
- \Uo\Rez 



(43) 



and 



G (z) = § exp 



C = ^|^exp(2|^|A m o) 



exp (2^' A 



pa y 



(44) 
(45) 
(46) 



It is necessary to carry out the complex time integrals. 



44 



The Zz integral. Here we have 



/ dzs exp | — ^ exp (— \Uq \ z) — \Uq \ Re , 
Jc 3 l D 



'C 3 

along the real axis of z, we can substitute 

u = exp (— \Uq \ z) 

and obtain 



LcM duexv ("t") = -WcJ dwexv( - w) 



This simply results in the constant 

D 



(47) 



The z 2 integral; The Z2 integral involves only the exponential of G (z 2 ) 
and the latter function does not contain a term proportional with Re z 2 
(comapre with (|43jl ). but only the exponential term. The function G (z 2 ) 
is very close to 1 as soon as 



U" V D 



Rez 2 > t ( s> ~ — In 



where is the Suzuki time for the region a 

Then, with very good approximation 

/ ef^ exp [G (z 2 )] = / dz 2 exp [G (z 2 ) 
ic 2 Jo 



t 
2 



The integral. This is similar to the integral over 23 and will imply again 



D 



as in (|4*7j) . 



Finally we note that only the constant Co relative to the intermediate 
point appears, and the one for the farthest point of the trajectory, a, is 
absent. 

There will be then a time integration over the moment of the first transi- 
tion, to? of an integrand that contains the time t arising from the integration 
on the contour C 2 , in the complex plane of the variable z 2 = t 2 — t\. This 
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variable represents the duration between the last transition to a and the first 
transition from a, i.e. is the duration of stay in a. This time integration will 
produce a term with t 2 . 

Replacing in the expressions resulting from the integrations on z\ and z% 
the constant Co from (|%H|) the factors x m and x p disappear. 

The result is 

*2> 







TJ"\ 




, 1/a (l^l 


^) V2 (K'l 


^') 1/2 


V2ttL) J 




2tt 


2tt 



x exp 



2^o -U a -U b 
D 



5 The probability 

The objective is to calculate P (b, |; 6, — |) on the basis of the above results. 
We have to sum over any number of independent pseudomolecules of the four 
species 

(bOb) < — ► -a b 
(60a) < — > a;, 
(a06) < — > a a 
(606) <— > -a a 

where the factor which is associated to each pseudomolecule is 

(K\uD 1/2 ( Uo - Ui \ . , 

a i — o ex P f; — , i = a,b 



2tt r \ D 

A particular path containing n (606) and m (aOa) has the form 

j.n+m+2p / \ TT"\ \ 1 / 2 

(_-] ) n + m n n +P n P+ m I 1 b 1 ) 

1 ' {n + m + 2p)\ b a \2ttD) 

It was taken into account that a path (6 — > 6) contains the same number p 
of (60a) and of (a06). 

5.0.1 Transition and self-transition probabilities at the equilib- 
rium points 

The results from the previous calculations are listed below. 

The probability that a particle initially at (6, — |) will be found 

at (6, §) 

t . t\ f\U"\\ 1/2 1 



P & > o' 6 ' ~o = o n 1 {°« + a ^ ex P H ( a « + a &)]} 

2 2 J \2-ku ) a a + o>b 
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The probability that a particle initially at (b, — |) will be found 

at (a, §) 



t . / N / A'" \ V2 



(^) ^7 ^ 1 " exp ["* ^ + 

\2ttVJ a a + a b 



The probability that a particle initially at (a, — |) will be found 

at (6,|) 

p (4 «■ 4) = ® V! ^ {1 - exp [ - ( <«• + a ' )l} 

The probability that a particle initially at (a, — |) will be found 

at (a, |) 

/ t t\ f\U"\\ l/2 1 
P^,-;a,--j = (^j __ { a 6 + a a expH(a a + «,)]} 

5.0.2 Transitions between arbitrary points 
The long time limit of the probability density 

The formula has been derived for the case where the two positions x and xq 
belong to the harmonic regions (a) and (b). This is because asymptotically 
these regions will be populated. 
Consider the case where 

x E (b) 
x G (a) 

The first trajectory contributing to the path integral is the direct connection 
between b and a. The contribution to K is 

redirect ( t M f S b0 + S 0a 

K x, -;x , -- ~ exp ' 



T u ' 1) r \ D 

Since the points of start and/or arrival x ,x are different of b, a, there is 
finite (i.e. non exponentially small) slope of the solution, and there is nomore 
degeneracy with respect to the translation of the center of the instantons. 
then there will be not a proportionality with time in the one-instanton term. 



47 



The terms connecting xq with x with onely one intermediate step in either 
b or a have comparable contributions to the action. 
Now the term with two intermediate steps 

Xq — > b — > a — > x 

In the very long time regime, can be factorized 



t 



t 



X 



/2tiD 



1/2 



if {1) (a,t';b,t") 



xKh arrn \ b,t ',xq, 



2ttD 



1/2 



This equation is independent of t' and t" as long as 



t K > | and f" 



t 



>^5 



Then 



if, 



harm 



•£) 2 ' a ' ^ 



V27T.D y 



1/2 



cxp 



t/(ar)-t/q 
2D 



There is also the approximate equality 



It results that all the contributions are included if in the expression of 

(x, |; Xo, — |) we replace the middle factor, which has become (a, |; fe, — |) 
by the full if (a, |; 6, — |). It is obtained 



t t 
P ( x, -; x , -- 



exp 



exp 



2L> 

t/(x)-t/g 

2D 



(rr// \ 1/2 
97h -^-{l-exp[-t(a a + a b )]} 
27rD / a a + a b 



The case where both x and x belong to the harmonic region (6) 

The same calculation shows that 



pL^Xo,- 1 - 



cxp 



U (x) - U b 
2D 



1/2 



i 



\2ttD J a a + a b 



{a a + a b exp [-t (a a + a b )]} 
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6 Results 



In the following we reproduce the graphs of the time- dependent probability 
distribution of a system governed by the basic Langevin equation. Various 
initialisations are considered, showing rapid redistribution of the density of 
presence of the system. The speed of redistribuition is, naturally, connected 
to the asymmetry of the potential. Each figure consists of a set of graphs: 

• the potential functions V(x), U(x) and W(x) 

• four probabilities of passage from and between the two equilibrium 
points a and 6; 

• the average value of the position of the system, as function of time, for 
the two most characteristic initializations 

The two parameters a and D take different values, for illustration. (We 
apologize for the quality of the figures. Better but larger PS version can be 
downloaded from |http : / /florin . spineanu . free . f r /sciarchive /topicalreview. ps | ) 
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Figure 5: Parameters: a = 0.4 and D = 0.001 
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Figure 6: Parameters: a = 0.4 and D = 0.005 
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Figure 7: Parameters: a = 0.4 and D = 0.05 
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Figure 8: Parameters: a = 0.6 and D = 0.005 
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Figure 9: Parameters: a = 0.6 and D = 0.05 
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Figure 10: Parameters: a = 0.75 and D = 0.005 
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Figure 11: Parameters: a = 0.75 and D = 0.05 
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X 

Figure 12: Parameters: a = 0.75 and D = 0.5 
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* 

Figure 13: Parameters: a = 0.9 and D = 0.05 



67 



Appendix A : Time variation for general form 
of the potential 

The stochastic motion is described in terms of the conditional probability 
P (x, t; Xi, 0) that the particle initially (t = 0) at Xi to be at the point x at 
time t. We will use the notation that suppress the as the initial time, 
P (x,t; Xi). The conditional probability obeys the following Fokker-Planck 
equation 

dP d rrTl . . „. d 2 P 

where the velocity function is here derived from the potential U 

F[x(t),t] = -U' (x) 

since there is no drive. The initial condition for the probability function is 

P (x, 0; x^ =6(x- x^ 

The simple form of the potential U (x) allows the introduction of two 
constants 

lo = |l/"(a; = 0)| 
wi = U" (x = 1) 

The following orderings are assumed 

D < uj 
D < ui 

The solution of the Fokker-Planck equation is given in terms of the fol- 
lowing path-integral 

exp [-U (x) /2D] T , , 
P (x, t] Xi) = P . L tt ) \'* K (x, t; Xi) 



where 



exp [-U{xi)/2D\ 



s_ 

D 



(A.l) 



where the functional integration is done over all trajectories that start at 
(xi,0) and end at (x,t). The action functional is given by 



S (x, t; Xi) = I dr 
'o 



1.2 

-x + W (x) 
4 
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with the notation 



W(x 



12 



In order to calculate explicitely the functional integral we look first for 
the paths x (t) that extremize the action. They are provided by the Euler- 
Lagrange equations, which reads 



SS_ 

with the boundary condition 



--x c + w (O = o 



x c (0) 
x c (t) 



■1 i 
X 



The first thing to do after finding the extremizing paths is to calculate 
the action functional along them, S c . After obtaining these extremum paths 
we have to consider the contribution to the functional integral of the paths 
situated in the neighborhood and this is done by expanding the action to 
second order. The argument of the expansion is the difference between a path 
from this neigborhood and the extremum path. The functional integration 
over these differences can be done since it is Gaussian and the result is 
expressed in terms of the determinant of the operator resulting from the 
second order expansion of S. Then the expression (IA.1)) becomes 



K = N- 



1 



det L (x c ) 



S c 

V2 ex P I ~d 



(A.2) 



where 



L (X r 



6 2 S 



5x (t)' 



1<P_ 
2d^ 



W" (X r 



(A.3) 



x(t)=x c (t) 

The constant is for normalization. 

In order to find the determinant, one has to solve the eigenvalue problem 



L (x c ) - A n 







(A.4) 



with boundary conditions for the differences between the trajectories in the 
neighborhood and the extremum trajectory 
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It can be shown that 



N 2 1 dE 



det L (x c ) 2nDx c (0) dx 



where E is the energy of the path x c (r). (see Dashen and Patrascioiu). 

There is a fundamental problem concerning the use of the equation (jA.2|) . 
It depends on the possibility that all the eigenvalues from the Eq. QA.4p are 
determined. However, a path connecting a point in the neigborhood of the 
maximum of U (x), 



is a kinklike solution. This solution (which will be replaced in the expression 
of the operator in Eq. (jA.4|) contains a parameter, the "center" of the kink, 
to- This is an arbitrary parameter since the moment of traversation is ar- 
bitrary. The equation for the eigenvalues has therefore a symmetry at time 
translations of to and has as a consequence the appearence of a zero eigen- 
value. Then the expression of the determinant would be infinite and the rate 
of transfer would vanish. Actually, the time translation invariance is treated 
by considering the arbitrary moment to as a new variable and performing a 
change of variable, from the set of functions <p n to the set (t , 4>n,n^o)- The 
trajectory that extremizes the action is a kinklike solution x (r, to) (is not 
exactly a kink since the shape of the potential is not that which produces 
the tanh solution) connecting the points : 



This treatment then consists of considering to as a collective coordinate (see 
Rajaraman and Coleman and Gervais & Sakita). The new form of 



where S is the action of the path x (r, to). The fact that the zero eigenvalue 
has been eliminated is indicated by the '. 




with a point in the neighborhood of one of the minima 




x(0,t ) = Xi 
x(t, to) = x 



Ea. (|A~2|) is 




(A.5) 
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The factor in the integral is approximately 

J2= 2^D [U " U {X)] 
Two limits of time are important. The first is 

2wo V D 
For the times 

< t < t a 

for a particle initially in the region of the unstable state, Xi ~ O (y/~D\ , the 

region around the stable state, ±1, or: x ~ ±1 +0 (v^D^ is insignificant for 

the calculation of the probability P (x, t; Xi). Then the expression (|A.2|) can 
be used. 

Much more important is the subsequent time regime 

t > t s 

where the particles leave the unstable point and the equilibrium state with 
density around the stable positions ±1 is approached. Then the general 
formula (I A. 51) should be used. 



The expression of the potential in general does not allow an explicit de- 
termination of the extremizing trajectories. Approximations are necessary. 

The energy E of the path x c (r) connecting the points (xi = 0, 0) and 
(x, t) can be approximated from the expression (see Caroli) 

dx' -pr - e (X, 0) 

where 

[* ( 1 

@o 0, y) = dx 



U' Q {x')\ \U'(x')\ 



The notations are 



Wo = \% 2 - jUS 

Uq (x) = —OJqX 

The last expressions are the harmonic approximation of Vq and the potential 
around the unstable point x — 0. 
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The additional term (x, 0) is the amount of time accounting for the 
deviation of V from the parabolic profile, the anharmonic deviation V — V . 
The result for this simpler case is 



u 1 

f [x, z; u, u) = 

where 



P (x, t; 0, 0) = ^^^VG exp (-G) 
Vtt \ U' (x)\ 



G = UJ ° %2 



D exp [2u T (x)] — 1 

T (x)=t + O Or, 0) 

Two time regimes are intersting. First, the very short time, less than the 
period of the particle in the potential, 

< t < 

the result reduces to the harmonic problem around x = 0. 
For longer times, 

^ <<t;St « 

the result is 

P (x, t; 0, 0) = -1=F' exp (-F 2 ) 

where 

F= 2 x exp [-cuqTo (x)} 

This is the scaling solution of Suzuki. 

The general expression of the function K in the case where we include the 
time regimes beyond the limits given above. It has the form of a convolution 

K= dt x (t ) K (x,ti,x m ) K (x m ,t ;x,i) 
Jo 



det L (x c ) 



1/2 



are obtained by the 



Here the approximations for the ratio N/ 

same method as before. The relation between the energy and the time is 
used for the two main regions: around the stable and the unstable (initial) 
positions 

t = dx'— —r - ©o (X m , Xi) 

[2(W + E )] 1/2 



X 



xm [2 {Wi + Ei)] / 
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